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Abstract: In this paper a finite difference/local discontinuous Galerkin method for the fractional 
diffusion-wave equation is presented and analyzed. We first propose a new finite difference method 
to approximate the time fractional derivatives, and give a semidiscrete scheme in time with the trunca¬ 
tion error 0((At) 2 ), where At is the time step size. Further we develop a fully discrete scheme for the 
fractional diffusion-wave equation, and prove that the method is unconditionally stable and convergent 
with order 0(h k+1 + (At) 2 ), where k is the degree of piecewise polynomial. Extensive numerical examples 
are carried out to confirm the theoretical convergence rates. 
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1 Introduction 


Fractional calculus, which might be considered as an extension of classical calculus, attracts 
much attention in recent decades. Fractional order partial differential equations (FPDEs) have 
been frequently used to solve many scientific problems in various fields, such as quantitative 
finance, engineering, biology, chemistry, hydrology, and so on 

However, analytical solutions for the majority of fractional partial differential equations, which 
are too complex and cannot expressed explicitly, are very difficult to be applied in the science 
and engineering, so it is a good choice to use numerical methods to finding numerical solutions 
for fractional partial differential equations, and has very important theoretical and practical 
significance. The existed methods solving the FPDEs include finite difference methods mum 
uni [131 [22U25U26U2ZI EZl |38l SSIIIZI S3] , finite element methods 0 HU El 1151 El E3 M , spectral 

methods [51 1251 [25], discontinuous Gakerkin methods [52152, homotopy perturbation method 
and the variational method [2J [29;, 32, [36] 153] [35. 55] . 

In this paper we consider the following fractional diffusion-wave equation 


d a u(x,t) 
dtH 

u(x, 0) = 


d 2 u(x,t ) 
dx 2 


= (x, t) £ [a, b] x [0, T], 


du(x, 0) 

Uo{ x ), — 7y t - =U i(x), x£[a,b\, 


( 1 . 1 ) 


where 1 < a < 2 is a parameter describing the order of the fractional time, f,uo,u\ are given 
smooth functions. We do not pay attention to boundary condition in this paper; hence the 
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solution is considered to be either periodic or compactly supported. 


The time fractional derivative in the equation (11.11) . uses the Caputo fractional partial deriva¬ 
tive of order a, defined as [TO] 



t > 0, 1 < ct < 2 


( 1 . 2 ) 


where T(-) is the Gamma function. 


The fractional diffusion-wave equation is obtained by replacing the first- or second-order time 
derivative of the classical diffusion or wave equation with a fractional derivative of order 1 < 
a < 2, and can be used to interpolate the diffusion equation and wave equation and model many 
of the mechanical responses and acoustics accurately. 


The rest of this paper is constructed as follows. In the section 2 some basic notations and the¬ 


oretic results are introduced. Then in section 3 we construct our finite difference/discontinuous 


Galerkin method for the fractional diffusion-wave equation, and stability and error analysis are 
given. Numerical results are presented in section 4, and the concluding remarks is included in 


the final section. 

2 Notations and auxiliary results 

In this section we introduce some notations and definitions that will be used later in the 
following sections. 

Let ft = [a, b] be a finite domain, and a partition is given by 


a = x i < X3 < ■■■ < x N ,i = b 
2 2 JV ' 2 


we denote the cell by L = \x ■ i, x il, for j = 1, • • • IV, and the cell lengths Ax,- = x-,i — 

J J 2 J J'2 

x- i, 1 < ? < N, h = max Ax,-. 

3 ~a — l<j<N J 

We denote by u + 1 and u~ 1 the values of u at x, +1 / 2 , from the right cell Ij + \ and from the 
left cell Ij, respectively. 

The piecewise-polynomial space is defined as the space of polynomials of the degree up to 
k in each cell Ij , i.e. 


V h k = {v:ve P k (Ij),x G Ij,j = 1, 2, ■ ■ ■ N}. 


For error estimates, we will be using two projections in one dimension [a,b], denoted by V, 
i.e., for each j, 



( 2 . 1 ) 


and special projection V ± , i.e., for each j, 



and 'P + o;(x + _ i) = cj(x-_ i) 
2 3 2 



and V w(x. i ) = uj(x, , i ). 



( 2 . 2 ) 
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For the above projections V and V^, we have [51E3 na HU 

||w e || + h\\i0 e \\ oo + h^\\u: e \\ Th < Ch k+1 , (2.3) 

where co e = Vco — lo or w e = — ui. 

The notations are used: the scalar inner product on L 2 (D) be denoted by (•,■)£>> and the 
associated norm by || • ||z>- If D = fl, we drop D. In the present paper we use C to denote a 
positive constant which may have a different value in each occurrence. 


3 The schemes 

In this section, we first present a finite difference method to approximate the time fractional 
derivatives, and then give the implicit fully discrete scheme with space discretized by the local 
discontinuous Galerkin method. Stability and convergence are detailed analysis. 


3.1 Time fractional derivative discretization 


We divide the interval [0, T] uniformly with a time step size At = T/M , M € IN, t n = 
nAt, n = 0,1, • • • , M be the mesh points. 

Let v(x,t) = 9u , and from the fact 

du(x,ti ) 3u(x,ti) - 4u(x,ti-!) + u(x,ti- 2 ) , „ 
v(x,U) = —^ = - — -+U. 

where the truncation error \r™\ < C(At) 2 , we can obtain 
d a u(x,t n ) 1 f tn dv(x,s) ds 


dt c 


1 

r (2 - 

a) 

i 


r (2 - 

a) 

i 


r (2 - 

a) 

(At) 2 

—a 

r(3- 

a) 

(At) 1 

—a 

r(3- 

a) 

(At) 1 

—a 

r(3- 

a) 

+ £ 

(' bn- 

i= 1 


bn— 

iv(. 


n —1 

___ / H 


dv(x,s) ds 


\cz —1 


»=o 

f ti+1 v(x,t i+ i) - v(x,ti) ds 

i =0 
n —1 

E" 




At (t n - s)' 

u(x,t i+ i) - v(x,ti) 


— + r 2 


+ P> 


i =0 


72—1 


- 6n_iu(x,t 0 )] + r% 


2At 


3u(x, tj) - 4tt(%, tj_i) + u(x, L_ 2 ) 
2A t 


where 


6o = l, bi = {i + iy~ a -i z ~ a ,i = 1,2,3, 


(3.1) 
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when i = 1, we take u(x, —1) = «(x,0) — Atu\(x) + C(At) 2 by Taylor expansion. 
Similar to the proof in [23],the truncation error \r%\ < C(At) 3 ~ a , so rj satisfied 

\r%\ < C(At) 3 ~ a . 


It is easy to check that 


bi > 0, i = 1 , 2 • ■ ■ , n. 

1 = bo > b\ > b 2 > • • • > b n , b n — > 0(n —> oo). 


(3.2) 


Substituting (13.11) into m, we have 


3 u(x,t n ) 


d 2 (x,t n ) 
dx 2 


n —1 

= ^(bn-i-i - b n _i)( 3(x, ti) - 4 u(x, ti- 1 ) + u(x, ti- 2 )) 

i =1 

+ 2Atb n -\v(x, t 0 ) + /3f(x , t n ) + 4 u(x, i) 

- u(x,t n - 2 ) + /3r%, 


where (3 = 2(At) Q T(3 — a). 

Let u k be the numerical approximation to u(x,tk), f n = f(x,t n ), the problem (11.11) can be 
discretized by the following scheme 


B 2 


n —1 


3u n - P-fop =^2( b n-i-i - K-i)(3u l - Au 1 1 + v l 2 ) 




+ 2A tb n -!V° + f3f n + 4 u n ~ l - u n ~ 2 , 


(3.3) 


where u 1 = u° — Atu±(x). We know 


|/3r£|<C(Af) 3 , 


however, by Taylor expansion we have 

| u(x, t- 1 ) — u(x, 0) + Atui(x)\ < C(At) 2 , 
therefore the truncation error is 0{At) 2 in scheme (13.31) . 


3.2 Fully discrete schemes 

In this subsection we present the fully discrete LDG scheme for the problem (11.11) based on 
the semidiscrete scheme (13.31) . 

We rewrite Eq. m as a first-order system: 

d a u(x 1 t ) .. . . 

P = u x , —— - p x = f(x,t). (3.4) 

Let G V k be the approximations of u(-,t n ),p(-,t n ), respectively, f n (x ) = f(x,t n ). We 

define a fully discrete local discontinuous Galerkin scheme as follows: find iiTp? £ V k , such 
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that for all test functions <f,w G Vjf, 


; [ ul<t)dx + f3( [ Vh^xdx - yZdPhfi 

Jn Jn 1=] 


) j+ i ~ C Ph<t> + )i-i )) 


n —1 


- 6„_j) / (3?4 - 1 + 2 )cj)dx + 2Atb n -i / 

Jn Jn 


(3.5) 


Z=1 

+ 4 f u!^~ l (f)dx — f u^~ 2 cj)dx + /3 [ f n (fdx, 

Jo, Jn Jn 

r r n 

/ plwdx + / u%w x dx - ^2((u%w -) j+ 1 - = 0, 

Jn Jn 2 2 

The initial conditions uf 1 , u° h ,v^ are taken as the L 2 projections of u(, — 1), u(, 0), u\ (, 0), respec¬ 
tively, 

/ uJ L 1 (f)dx= / Vu(x,—l)(j)dx= / u(x, — l)</>dx, 

Jn 

/ Ufrtfidx = / Vu(x, 0)<pdx = / uo{x)(j)dx , (3.6) 

Jo Jo Jo 

/ Vfrjidx = / Vui(x, 0)cf>dx = / ui(x)<j)dx, Vu G Vjf . 

The “hat” terms in (13.51) in the cell boundary terms from integration by parts are the so- 
called “numerical fluxes”, which are single valued functions defined on the edges and should be 
designed based on different guiding principles for different PDEs to ensure stability. It turns 
out that we can take the simple choices such that 


K = W) , Ph = ( Ph) + ■ 


(3.7) 


We remark that the choice for the fluxes (|3.7I) is not unique. In fact the crucial part is taking 
and p^ from opposite sides |50j 15] . 


3.3 Stability and Convergence 

In order to simplify the notations and without lose of generality, we consider the case / = 0 
in its numerical analysis. 

Theorem 3.1. For periodic or compactly supported boundary conditions, the fully-discrete 
LDG scheme is unconditionally stable, and there exists a positive constant C depending 

on u , T, a, such that 


IKII < C(\\u ° h || + At|Mx)||), n = 1, 2 • • • , M. (3.8) 

Proof. Taking <f = u^,w = (Jp 1 ^ in scheme (13.51) . we obtain 

N 

3IKII 2 + «ll 2 + j®, + > - 'I'K.Kyi + e«, 

1=1 
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(3.9) 


71—1 


y2(b n -i -1 - b n -i) / (3 u\ - 4 u\ 1 + 2 )<ete + 2Af6„_i / v^u^dx 
■ _ 


2=1 


+ 4 [ nl~ x uldx - [ u n h ~ 2 u n h dx, 
J Q, J 0, 


where 


=K)-«)“-KKr-<W' 

©K,K) =(pST «)" - (Ph) + K) + -®«h)“ +^K) + - <{pir 

+ K(Ph) + - 

If we take the fluxes (13.71) . after some manual calculation, we can easily obtain @(u^,p^) = 0. 
Then based on the equation (13.91) . we can get 


71—1 


77. 11 2 


+ P\\Phf =^2{b n -i-l ~b n -i 

i= 1 


) [ ( 3 < 

Jn 


- 4 < 1 + < 2 K d ® 


+ 2Atb n -i / v^ii^dx + 4 v,fr l u r J L dx 
Jn Jn 

- [ ul~ 2 uldx 
Jn 


n— 1 


< 


^2(p n -i -1 - & n-i)(3|K|| +4||< X || + ||u 


2 —211\ iL.n | 

l h 


2—1 


+ 2 At 6 B _ 1 KHK||+ 4 K- 1 ||K|| 

+ K“ 2 IIKII, 


that is 


71 — 1 


< 


J2( b n-i-l - «>n-i)(3|K|| +4||< 1 || + |K 


2-2 I 


2—1 


+ 2At6 n _i|K|| + 4 IK 


n—1 1 


(3.10) 


+ IK“ 2 I 


We will prove the Theorem 13.II bv mathematical induction. When n = 1, we can obtain 

3|K|| < 2A2||ug|| +4|K|| + IK 1 !! (3-11) 

Notice that 

[ u^vdx = f P(«(x,0) - At Ul (x))vdx = j u° h vdx-At j u x (x)vdx, 

I A %J I A _/" A *J J A 


for any v G Vfi. Taking v = u h 1 , we can obtain 


u h 11/ = / uj l u h 1 dx — At / u\(x)u h dx 
3 A Jij 

^ ll^hll/,- + ~^\\ u h L \\‘ij + (^) 2 ll' u i( x )ll/j + jll u / i 1 || 2 i ) 


6 


summing over j from 1 to N, we can get 


IK 1 !! A c '(ll“°ll + A *IM®)ID- ( 3 - 12 ) 

Similar to the proof of (13.121) . we can easily obtain 

Kll < HxiWII- (3.13) 

By using (I3.11IUI3.12D and (13.131) . it is easily to know that there exists a positive constant C, 
such that 


Kll<C(K||+At|Mz)||). (3.14) 

Now suppose the following inequality holds 

IKII < C{\\u° h \\ + At|| Ul (x)||), m = 2,3 • • • K, (3.15) 

we need to prove ||n^ +1 || < (7(||u°|| + At ||rti(x)||). 

Let n = K + 1 in the inequality (13.101) . we can obtain 

K 

3||»f +1 ll <£( 6 K-<- i «+i-<K 3 K!l + 4 K-'II + IK -2 II) 

i=1 

+ 2AtbK\\vh\\ + 4||n^ r || 

+ IK" 1 H- 

Using ()3.121) . (13.131) and (13.151) . we can obtain the following inequality easily 

||wf +1 || < C{\\u° h \\ + At||ui(x)||). 

This finishes the proof of the stability result. □ 

Theorem 3.2. Let u(x,t n ) be the exact solution of problem which is sufficiently smooth 

such that u £ H m+1 with 0 < m < k + 1. Let u^ be the numerical solution of the fully discrete 
LDG scheme m>, then there holds the following error estimate: 

|| u(x, t n ) — u'f || < C(h k+1 + (At) 2 ), n = 1, ■ ■ • , M, (3.16) 

where C is a constant depending on u,T,a. 

Proof. By Taylor expansion we know 

| u(x, t- 1 ) — u(x, 0) + Atui(x)\ < C(At) 2 , 

here C is a positive constant depending on u. Then by using the property (12.31) . we can obtain 
the following estimate which will be used later, 

\\u(x,t-i) — uf l \\ < C((At) 2 + h k+1 ). (3-17) 
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It is easy to verify that the exact solution of PDE (11.11) satisfies 

I f \'' 

3 / u(x,t n )cj)dx + /3( p(x,t n )(/> x dx - f n )0“) -,i - {p{x ,t n )4> + ) ■ 1 )) 

Jn Jn “( 3+2 32 

n— 1 r. 

- b n _i) / (3u(x,ti) - 4u(x,ti-i) + u(x,ti- 2 ))(j)dx 

Jn 

/ v(x, to)(j)dx + 4 / u(x,t n -\)4>dx 

J Q J n 

h t n ~ 2 )(j)dx + /? / f(x,t n )(/)dx + /3 / r%4>dx, 

Jn Jn 


n —1 

£< 

+ 2At6 n _i 


(3.18) 


- / U{X,t n -<2 

Jn 


f f J' 

/ p(x,t n )wdx + / u(x,t n )w x dx - y^((u(x,t n )w~) ■, i - (u(x, f n )u; + ) ,_i) = 0, 

In -In 3 2 3 2 

Vv,r] G H 3 (Ij), for j = 1, ■ ■ ■ JV. 

Denote 

e" = u(x,t n ) -u% = V~e ™ - (V~u(x,t n ) - u(x,t n )), 
e p = p(®, - Pft = V + e™ - (V + p(x, t n ) - p(x, t n )). 

Subtracting (13.51) from (13.181) . and with the fluxes (13.71) we can obtain the error equation: 

r r 33 

3 e^dx + fK e;^dx-Y.(((e;) + r) j+ L-(K) + <l> + ) j . i )) 

Jn Jn • 1 2 2 


(3.19) 


71—1 


^2(bn-i -1 - 6n-i) / (3< - 4e^ 1 + 4 2 )</Wx 


2=1 

-4 


/ e" Vdx 
Jn 

+ / e^~ 2 (j)dx + (3 r^cfrdx + / e™wdx + / e™w x dx 

J n J n J n J n 


(3.20) 


AT 


- - ((^r^-i) = o. 

1=1 

Using (|3.19l) . the error equation (13.201) can be written as follows: 


) i+ j - ((r + 4) + ^ + ) i -i)) 


I [ V e™4>ix + P( ( P+eJ0 l( ii - V(((P + e”) + <^ 

Jn Jn . =1 

r r n 

+ / P+e>ds+ / iP- e > 3: dx-^(((iP- e -)- U ;-) J . + i - ((iP-e^)-^)-_i) 
in in “[ 2 2 


71—1 


J=1 

- 6n-<) / (3PX - 4p-e(r 1 + V-(J- 2 )<j)dx 

i= l ^ 

+ 4 f V~e™- 1 4>dx — f V~e™~ 2 <j)dx — j3 f r’Jipdx 
J n J n J n 
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(3.21) 


71 — 1 


- y~](frn-i-i - K-i) / (3(p u(x,ti) - u(x,ti)) - 4(P u(x, ti- 1 ) - tt(x, t j_i)) 
i=l 7n 

+ (V~u(x,ti- 2 ) - u{x,ti-2)))4>dx - 4 / (P~u(x,t n -i) - u(x,t n -i))(/)dx 

Jn 


+ (V u(x,t n - 2 ) - u(x,t n - 2 ))(j}dx 
Jn 


+ 3 (V u(x,t n )-u(x,t n ))(j)dx + P( (V + p(x,t n ) -p(x,t n ))(j> x dx 
O JO 

N 

~ ^(((CP + p(Mra) -P(^,tn))) + 0") j+ i - (({V + p{x,t n ) -p(x,^ n ))) + ^ + ) j _l)) 
i=i 

+ / CP + pOMn) -p(x,t n ))wdx + / (V~u(x,t n ) — u(x,t n ))w x dx 
JO JO 

N 

- y^({(V~u(x,t n ) - u(x,t n ))~w ~) j+ 1 - ((V~u(x,t n ) - u(x,t n ))~W + ) j _l). 

3 = 1 

Taking the test functions (j> = V~e™,w = j3V + ep in (13.211) . using the properties (l2.lD - d2.2f) . 
then the following equality holds, 

3 f (V~e™) 2 dx + P [ ( V + ep) 2 dx 
Jn Jn 

n —1 „ 

= V(b n _i_i - ifc-i) / ( 3 P-ej, - diP-e 4 " 1 + V~j- 2 )V-e n u dx 

i =i 


■f 

Jn 


+ 4 / V-e^PS^dx - I V~eVT*P~eVidx 


"U ' U 


I 

Jn 


->—„n—2 -n- 


"22 u 


71—1 


- & n _j) / (3(7 ? - u(x,ti)) - 4(P n(x,ti_i) 


2 — 1 


- n(x,ti_i)) + ("P u(x,ti- 2 ) - u(x,ti- 2 )))V e^dx - 4 

[ (V~u(x,t n -i) - u(x,t n -i))V~e™dx - /3 f r^V^e^dx 
Jn Jn 

+ / (P~u(x,t n - 2 ) - u(x,t n - 2 ))V~e™dx 
Jn 

+ 3 (V~u(x,t n ) - u(x,t n ))V~e™dx 
Jn 

+ 13 [ (V + p(x,t n ) - p(x,t n ))V + epdx. 


Therefore, we obtain 


71—1 


-2—2 I 


ll^-OII <2^(6n-i-r -6n-i)(3r-clll +4||P-eni + H^e 

2—1 

+ 4||p-eT 1 || + r-er 2 ||+/3||r 3 "|| 

72—1 

+ y^(b n _^_i - 6 n _i)(3||T , ”n(x,t i ) - rt(x,^)|| 

2 — 1 

+ 4||T ,_ n(x,t i _i) - n(a:,ti_i)|| + ||T ,_ n(a;, U- 2 ) - u(x,U - 2 )||) 
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(3.22) 


+ \[P u(x,t n - 2 )-u(x,t n - 2 ) ||+4||P u(x,t n -i) - u(x,t n -i)\\ 

+ 3|| V~u(x,t n ) - u(x,t n )II + y/J}\\V + p(x,t n ) -p(x,t n )||. 

We prove the error estimates (13.161) by mathematical induction. First, we consider the case 
when n = 1. From (13.201) and (|3.23l) , we know 

WV-eDW <4|p-e°|| + || e - 1 ||+/3||r 3 1 || 

+ 4\\V~u(x,to) - u(x,t 0 )\\ (3.23) 

+ 3\\V~u(x, ti) - u(x,ti)\\ + y/]3\\V + p(x,ti) -p(x,t 1 )\\ 

Notice the facts that 

V~e °=0, He" 1 !! <C(h k+1 + (At) 2 ), ||r£|| < C{Atf~ a , 

and the property (12.31) . we can obtain 

||P“ei|| < C(h k+1 + (At) 2 ). (3.24) 

Next we suppose the following inequality holds 

\\V~e™\\ < C{h k+l + (At) 2 ), m = 1, 2, • • • K. (3.25) 

When n = K + 1, from the equation (13.231) . we can obtain 
K 

l|Wef +1 )|| <]>>*_, - 6 K+ i_i)(3||T>-<|| + 4||p-eJ,- 1 || + ||P^r 2 ||) 

i— 1 

+ M\V-eZ\\ + \\V-eZ- 1 \\+P\\r« +1 \\ 

K 

+ ^ ~2(b K -i - b K +i-i){3\\V~u(x, ti) - u(x,ti) || + A\\P~u{x, t*—i) - u(x,U- 1 )|| 

2— 1 

+ || V~u(x,ti- 2 ) - u(x,tj_ 2 )||) + 4||P _ u(x,tft:) - u(x,t^)|| 

+ \\V~u(x,t K -i) - u(x,t K - 1 )|| 

+ 3||7 ?_ u(a:,tii: + i) -u(x,t K + 1 )|| + V^II^’pOM-K'+i) — p(£c, t^+i)|| 

K 

< Y,{b K -i ~ b K +i~i)C(h k+1 + (At) 2 ) 

2=1 

+ 4C(h k+1 + (At) 2 ) + C(h k+1 + (At) 2 ) + C(At) 3 

K 

+ - 6^ + i_i)C/i fc+1 + Ch k+1 + C(At)f h fc+1 . 

2=1 

Similar to the proof of (13.241) . we can obtain the following result immediately 

\\V~e^ +1 \\ < C(h k+1 + (At) 2 ). 

Thus Theorem 13.21 follows by the triangle inequality and the interpolation property (12.31) . □ 
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4 Numerical examples 


In this section, we present numerical experiments of the presented finite difference/local dis¬ 
continuous Galerkin method to the fractional diffusion-wave equation to verify the error esti¬ 
mates in Section 3. 

Example 4.1. Consider the following fractional diffusion-wave equation 


d a u(x,t) 
dt “ 

u(x, 0) = 0, 


d 2 u(x 1 t) 
dx 2 

du(x, 0) 
dt 


f(x,t), (x, t) G [0,1] x [0,1], 

= 0, x € [0,1]. 


(4.1) 


Choose a suitable right-hand-side function f(x,t ) such that the exact solution is 


u(x, t ) = t 2 sin(27rx) 


. Table 1-4 display the errors in L 2 -norm and L°°-norm at T = 1 and convergence orders 
in space for piecewise P k polynomials for several values of a : 1.2,1.4,1.6,1.8, with time step 
At = 1/1000. Obviously the (k + l)-th order of accuracy in space are observed, which is in 
agreement with the theoretic results. 

In order to investigate the temporal accuracy of the proposed method, we fix the space step 
h = 1/200. Table 5 show that the errors in L 2 -norm and id-norm attain the second-order 
convergence in time. The results are consistent with our theoretical results in Theorem 13.21 

In Figure Q] we plot the the approximate solution of the three order on the uniform mesh with 
100 cells and the exact solution at T = 1 to show the performance of the presented scheme. We 
can see that the method is very effective and is a good tool to solve such problems. 


5 Conclusion 

In this work, we have presented a finite difference/local discontinuous Galerkin method for the 
fractional diffusion-wave equation. We first propose a finite difference method to approximate 
the time fractional derivatives when 1 < a < 2, and then give a fully discrete scheme and prove 
that the scheme is unconditionally stable and convergent. The extensive numerical example 
fully confirm the theoretic analysis. 
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Table 1: Spatial accuracy test for the time-fractional diffusion-wave equation (|4.1D using 
piecewise P k polynomials, a = 1.2, A t = gggg,T = 1. 



N 

L 2 -error 

order 

L°°-error 

order 


5 

0.265109983989909 

- 

0.623532065154133 

- 


10 

0.129308265170869 

1.03 

0.313595441923762 

0.98 

pO 

20 

6.425848539792525E-002 

1.01 

0.157010259108059 

1.00 


40 

3.208010396964848E-002 

1.00 

7.853125916762804E-002 

1.00 


80 

1.603391954016398E-002 

1.00 

3.926891940062052E-002 

1.00 


5 

6.736979744152280E-002 

- 

0.249880240379731 

- 


10 

1.695495641564284E-002 

1.99 

6.468476942047330E-002 

1.95 

P 1 

20 

4.245128618901225E-003 

2.00 

1.631233665625631E-002 

1.99 


40 

1.061671545606701E-003 

2.00 

4.103801542149732E-003 

1.99 


80 

2.654420510184321E-004 

2.00 

1.027598756713433E-003 

2.00 


5 

6.682959934132981E-003 

- 

3.174066978350254E-002 

- 


10 

8.506996364720942E-004 

2.97 

3.971254358826398E-003 

3.00 

p2 

20 

1.068204883018595E-004 

2.99 

5.116352220441455E-004 

2.96 


40 

1.336779544959365E-005 

3.00 

6.443554411463790E-005 

2.99 


80 

1.672333408099435E-006 

3.00 

8.069512789853388E-006 

3.00 


Table 2: Spatial accuracy test forthe time-fractional diffusion-wave equation EH using 
piecewise P k polynomials, a = 1.4, At = , T = 1. 



N 

L 2 -error 

order 

L°°-error 

order 


5 

0.265002133818103 

- 

0.623263433950011 

- 


10 

0.129296240434182 

1.03 

0.313564577675057 

0.98 

pO 

20 

6.425702562008787E-002 

1.01 

0.157006484631294 

1.00 


40 

3.207992382929516E-002 

1.00 

7.853079254830535E-002 

1.00 


80 

1.603389759999361E-002 

1.00 

3.926886254925913E-002 

1.00 


5 

6.736429618528465E-002 

- 

0.249848519443165 

- 


10 

1.695467033281010E-002 

1.99 

6.468306250704692E-002 

1.95 

P 1 

20 

4.245111564609582E-003 

2.00 

1.631223051257358E-002 

1.99 


40 

1.061670492069727E-003 

2.00 

4.103794461966181E-003 

1.99 


80 

2.654419852360282E-004 

2.00 

1.027597816101400E-003 

2.00 


5 

6.682553306760190E-003 

- 

3.173870305522863E-002 

- 


10 

8.506873719172307E-004 

2.97 

3.971185742546351E-003 

3.00 

p2 

20 

1.068201082305646E-004 

2.99 

5.116330685330400E-004 

2.96 


40 

1.336778230870020E-005 

3.00 

6.443547673697195E-005 

2.99 


80 

1.672322693664752E-006 

3.00 

8.069510669296185E-006 

3.00 
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Table 3: Spatial accuracy test forthe time-fractional diffusion-wave equation (14.11) using 
piecewise P k polynomials, a = 1.6, At = = 1. 



N 

L 2 -error 

order 

L°°-error 

order 


5 

0.264972688983567 

- 

0.623189941325662 

- 


10 

0.129291634839670 

1.03 

0.313552730827037 

0.98 

pO 

20 

6.425644681449535E-002 

1.01 

0.157004984816329 

1.00 


40 

3.207985524403432E-002 

1.00 

7.853061451680375E-002 

1.00 


80 

1.603389101024813E-002 

1.00 

3.926884544112283E-002 

1.00 


5 

6.736317907360344E-002 

- 

0.249838735643862 

- 


10 

1.695461845352325E-002 

1.99 

6.468254410369267E-002 

1.95 

P 1 

20 

4.245108546494545E-003 

2.00 

1.631218386012478E-002 

1.99 


40 

1.061670300965039E-003 

2.00 

4.103777047141932E-003 

1.99 


80 

2.654419703069185E-004 

2.00 

1.027582241539760E-003 

2.00 


5 

6.682478515070153E-003 

- 

3.173831828892504E-002 

- 


10 

8.506852141270998E-004 

2.97 

3.971173620498153E-003 

3.00 

p2 

20 

1.068200376865845E-004 

2.99 

5.116326863045966E-004 

2.96 


40 

1.336774817234452E-005 

3.00 

6.443546414282400E-005 

2.99 


80 

1.672066664180233E-006 

3.00 

8.069510170348080E-006 

3.00 


Table 4: Spatial accuracy test for the time-fractional diffusion-wave equation m using 
piecewise P k polynomials, a = 1.8, At = iggo,T = 1- 



N 

T 2 -error 

order 

L°°-error 

order 


5 

0.265395156138938 

- 

0.624238212597400 

- 


10 

0.129328509614864 

1.03 

0.313647188207608 

0.98 

pO 

20 

6.426055758264328E-002 

1.01 

0.157015597179358 

1.00 


40 

3.208037034587025E-002 

1.00 

7.853194655385835E-002 

1.00 


80 

1.603396389418712E-002 

1.00 

3.926903380844882E-002 

1.00 


5 

6.736868534645306E-002 

- 

0.249898961167015 

- 


10 

1.695486972560303E-002 

1.99 

6.468562726677418E-002 

1.95 

P 1 

20 

4.245123063266223E-003 

2.00 

1.631242414011957E-002 

1.99 


40 

1.061671210932901E-003 

2.00 

4.103841431837951E-003 

1.99 


80 

2.654420391246725E-004 

2.00 

1.027635651727810E-003 

2.00 


5 

6.682807553983295E-003 

- 

3.174004623627680E-002 

- 


10 

8.506954178648776E-004 

2.97 

3.971230808741821E-003 

3.00 

p2 

20 

1.068203742828856E-004 

2.99 

5.116345075529710E-004 

2.96 


40 

1.336790404833714E-005 

3.00 

6.443552350912754E-005 

2.99 


80 

1.673232029466757E-006 

3.00 

8.070278651961527E-006 

3.00 
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Table 5: Temporal accuracy test for the time-fractional diffusion-wave equation m using 
piecewise P 2 polynomials. N = 200 



At 

L 2 -error 

order 

L 1 -error 

order 


0.05 

2.315895458864733E-006 

- 

2.078028449959208E-006 

- 


0.04 

1.452896981500066E-006 

2.09 

1.301714631618507E-006 

2.10 

a = 1.1 

0.03 

8.051939547113251E-007 

2.05 

7.198086134684791E-007 

2.06 


0.02 

3.394311162501219E-007 

2.13 

3.023980936139506E-007 

2.14 


0.05 

1.550654829179151E-004 

- 

1.396164504933172E-004 

- 


0.04 

9.400505980906901E-005 

2.24 

8.464278635512412E-005 

2.24 

a = 1.8 

0.03 

5.271466986243190E-005 

2.01 

4.746819651138549E-005 

2.01 


0.02 

2.314548304581403E-005 

2.03 

2.101160952907237E-005 

2.01 


exact solution when T=1 


x 



0.8 


0.6 


0.4 


0.2 


0.0 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.9 


a 



1.000 

0.7500 

0.5000 

0.2500 

0 

-0.2500 

-0.5000 

-0.7500 

-1.000 


numerical solution when T=1 


x 




1.010 

0.7575 

0.5050 

0.2525 

-1.110E-16 

-0.2525 

-0.5050 

-0.7575 

-1.010 


1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

a 


Figure 1: The exact solution is contrasted against approximate solution obtained on the 
uniform mesh with 100 cells using P 2 elements for Eo. dl.ll) when T=l. 
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